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ABSTRACT 

We have conducted VI CCD photometry of the two open clusters NGC 
1245 and NGC 2506 using the CFH12K CCD camera. Our photometry 
covers a sky area of 84' x 82' and 42' x 81' for the two clusters, respectively, 
and reaches down to V ~ 23. We derived the physical parameters using 
detailed theoretical isochrone fittings using \ 2 minimization. The derived 
cluster parameters are E{B-V) = 0.24±0.05 and 0.03±0.04, (V-M v ) = 
12.25 ±0.12 and 12.47±0.08, age(Gyr) = 1.08 ±0.09 and 2.31 ±0.16, and 
[Fe/H]= -0.08±0.06 and -0.24±0.06, respectively for NGC 1245 and NGC 
2506, We present the luminosity functions (LFs) of the two clusters, which 
reach down to My 10, and derive mass functions (MFs) with slopes 
of T = -1.29 for NGC 1245 and T = -1.26 for NGC 2506. The slopes 
are slightly shallower than that of the solar neighbourhood, implying the 
existence of dynamical evolution that drives the evaporation of the low- 
mass stars in the clusters. 

Key words: (Galaxy:) open clusters and associations: individual: - meth- 
ods: observational - techniques: photometric 
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1 INTRODUCTION 

Open clusters are located in the Galactic disc as Popu- 
lation I, and it is important to investigate their prop- 
erties and the spatial distribution of the Galaxy to 
understand the Galaxy structure and evolution, espe- 
cially the formation and dynamical evolution of the 
Galactic disc. Because stars in the Galactic disc are 
originated in evaporated stars from open clusters, the 
luminosity function of open clusters is an important 
parameter for understanding their formation and evo- 
lution. We can easily obtain the LF from CCD obser- 
vations, but an analysis of the bF is not easy because 
the observed bF of open clusters reflects the dynam- 
ical evolution which depends on the cluster's age as 
well as the initial mass function (IMF). 
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It is necessary to have a deep photometry over the 
entire cluster field to understand the FF of old open 
clusters down to the magnitudes of the evaporating 
stars because they are much fainter than the turn- 
off stars and are supposed to be located around the 
tidal radius of the clusters. But there are only a few 
open clusters of which photometry is deep and wide 
enough to study the effect of dynamical evolution on 
the structure of open clusters. Thus, we selected open 
clusters which are old enough for the low mass stars 
to evaporate from the clusters by dynamical evolution 
and rich enough suitable for a detailed analysis of the 
spatial distribution of low mass stars. We also took 
into account the distance of open clusters to ensure 
that our photometry reaches up to Mv ~ 10. 

The old open cluster NGC 1245 is located in the 
Perseus arm at a Galactic latitude of —8°. 93. Us- 
ing photo graphic photometry data from th e Schmidt 
telescope. iPeruanskii fc Rvadchenkd (1 19781 ') discussed 
the radial density distribution of stars with a limiting 
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magnitude of B = 17 and a radius of r « 83'. They 
presented cluster boundaries with a nucleus (r < 13') 
and corona (13' < r < 6 1'). ICarraro fc Patatl i| 19941 ) 
and ISubramaniaml (|2003l ) found a s hallower slope o f 
the mass function of NGC 1245 than ISalpetei (1 19551) . 
Using BV CCD photometry, ISubramaniam |2003) 
presented the surface number density with a limit- 
ing magnitude of V = 17 and a radius of r ~ 6'. 5. 
They also indicated a lac k of stars in the cluster 
centre. iBurke et al. I l|2004h derived physical param- 
eters and a core radius of R c — 3'.1 from BVI CCD 
photometry. They determined an optical absorption 
of Ay = 0.68 magnitude and a distance modulus 
(V — My)o = 12.27; they also reported that the 
cluster has no differential reddening, in contrast to 



previous studies s uch as ICarraro fc Patatl l|l994h and 
IWee fc Led i| 19961 ). 

The old open cluster NGC 2506 also has rich 
members located at a Galactic latitude of +9°. 94. A 
radial density profi le of the proper mo t ion m embership 
stars obtained bv lChiu fc van Altenal (|l98 ll ) indicated 
that bright stars ar e more concen t rated than faint 
stars in this cluster. iMarconi et alj (|l997i ) suggested 
that more than 20% of the main s equence stars are 
binary . From a photometric study bv lHenderson et al.l 
(2007), the distance of this cluster is known to be 3.4 
kpc and its age is 1.79 Gyr. 

This paper presents the deep VI CCD photom- 
etry of the two old open clusters NGC 1245 and 
NGC 2506, covering the entire cluster fields including 
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the surrounding field regions for field star correction, 
along with the derivation of the cluster parameters 
such as age, metallicity, distance, and interstellar red- 
dening. We will discuss the dynamical structures of 
the open clusters in a forthcoming paper (Lee, Kang 
& Ann, in preparation), but here we present the phys- 
ical properties, luminosity functions, and mass func- 
tions of the two open clusters along with their physical 
properties. 

The present paper is organized as follows. In §2, 
we describe the observations and data reductions. We 
present the physical parameters of the clusters in §3, 
and describe the luminosity functions and mass func- 
tions in §4 and §5, respectively. A summary is given 
in the final section. 



2 OBSERVATION 

We selected target clusters best suited for studying the 
dynamical properties of open clusters which is the sub- 
ject of the second paper (?). The two target clusters 
NGC 1245 and NGC 2506 are located at a relatively 
high Galactic latitude to reduce field star contamina- 
tion and are old enough for the evaporation of low 
mass stars to be prevalent due to dynamical evolu- 
tion. Also, they are at proper distances for optimum 
sky coverage and the deep photometry to reach up to 
My ~ 10. The basic known properties of the selected 
clusters are listed in Table 1. 




2.1 Observation 

We conducted deep VI CCD photometry using 
the CFH12K mosaic CCD camera mounted at the 
prime focus of the Canada-France-Hawaii Telescope 
(CFHT). The effective aperture of the CFHT is 3.6 
m, and the CFH12K mosaic CCD camera consists of 
12 CCDs of 2048 x 4096 pixels. The CFH12K field of 
view is 42' x 28' with pixel size of 0."206. 

For better sky coverage, we observed 6 regions for 
NGC 1245 and 3 regions for NGC 2506. The field of 
view of the combined images was 84' x 82' for NGC 
1245 and 42' x 81' for NGC 2506 (see Fig. 1 and Fig. 
2). The observations were performed on the night of 
December 6, 2001, at CFHT. The observation log is 
given in Table 2. 



2.2 Data Reduction 

We followed the standard data reduction proce- 
dures, including bias subtraction, overscan correction, 
and trimming and flat-fielding using IRAF/CCDRED 
(|Stetson|[l987h . We used the twilight sky flat fields 
obtained before and after the observations for flat- 
fielding. In order to increase the signal-to-noise ra- 
tio and to reach low limiting magnitudes, we took 
multiple exposures of each region. These images were 
aligned and median-combined into a master image for 




Figure 2. Observed fields (42' X 81') of NGC 2506. We 
combined 3 V-filter images, where one CCD image has a 
size of 42' X 28', to make the mosaic image. North is up 
and east is to the left. 



each region using images of both the long and short ex- 
posures. For the photometry calibration, we obtained 
calibration images for each chip of the CFH12K. We 
used 20-40 isolated bright stars to derive the point- 
spread function (PSF) of each c alibration imag e, and 
we performed PSF photometry (|Stetsonlll987h using 
IRAF/DAOPHOT to obtain the instrumental magni- 
tudes. 

For standard star transfor mation in eac h band, 
we observed the SA 98 region (|LandoltH l992) on the 
same night. The standard star observations were care- 
fully designed to ensure that the maximum number of 
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Table 1. Basic data of the selected clusters, extracted from the WEBDA web site 
(http:/ /obswww.unige.ch/webda/ ) 



Cluster 


a 


S 


I 


b 


Distance 


E(B - V) 


log(age) 


name 


J2000.0 


J2000.0 






(kpc) 






NGC 1245 


03 h 14 m 42 s 


+47 14 12 


146° .65 


-8°. 93 


2.88 


0.300 


8.704 


NGC 2506 


08 ,t 00 m 01 s 


-10 46 12 


230°. 56 


+9°. 94 


3.46 


0.081 


9.045 



Table 2 . Log of observations 



Cluster 


Region 


V Filter 




/ Filter 




Seeing 


NGC 1245 


North-East 


140 s x 4, 5 s x 


4 


90 s x 4, 5 s x 


4 


l."6 




East 


140 s x 4, 5 s x 


4 


90 s x 4, 5 s x 


4 


l."2 




South-East 


140 s x 4, 5 s x 


4 


90 s x 4, 5 s x 


4 


l."l 




North- West 


140 s x 4, 5 s x 


1 


90 s x 4, 5 s x 


4 


l."3 




West 


140 s x 4, 5 s x 


1 


90 s x 4, 5 s x 


4 


l."5 




South- West 


140 s x 4, 5 s x 


1 


90 s x 4, 5 s x 


4 


l."3 


NGC 2506 


North 


140 s x 3, 5 s x 


3 


90 s x 4, 5 s x 


3 


l."2 




Centre 


140 s x 3, 5 s x 


3 


90 s x 4, 5 s x 


3 


l."0 




South 


140 s x 3 




90 s x 4 




l."2 



0.2 

> 



< 

-0.2 

i 

:> o 

< -0.2 
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Figure 3. Photometric residuals of the transformation 
to the standar d system. We used standard stars in SA 98 
iLandoltlll992T) . Standard deviations of the VI photometry 
are &av = 0.018 and cta(v-I) = 0.028. 
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Figure 4. Incompleteness correction factor A as a function 
of V magnitude. 



stars were imaged in each chip of the CFH12K mosaic 
CCD. The number of observed standard stars was 8- 
19 for each chip. We derived the standard magnitudes 
of the programme stars by using the following equa- 
tions given on the CFHT web page: 

V = v - 0.12(X - 1) + 0.0U(V - I) + z v (1) 
I = i - 0.04(X - 1) + 0.050(1/ -I)+ Zi (2) 

where V and I are the standard magnitudes and v 
and i are the instrumental magnitudes. The terms z v 
and Zi are the zero points of the photometry that were 
determined from the photometry of the standard stars 
in SA 98, and X is the airmass. Fig. 3 shows the resid- 
uals of the standard star transformation in the V and 
/ bands. 

To correct the incompleteness of the photometry, 
we generated artificial stars on the CCD frames us- 
ing the addstar task in the IRAF/DAOPHOT utility 
and performed reduction in the same manner as the 
observed frames. The additional artificial stars were 
randomly generated in each unit-magnitude bin on all 
original images. To prevent crowding in the images, 
the number of artificial stars were controlled to be 
less than 10% of the detected stars, and we repeated 
this for 40 trials. The incompleteness factor (A) of the 
photometric data was determined in the central and 
outer regions of each cluster. Finally, the incomplete- 
ness correction factor was obtained from A = N r /N a , 
where N a is the number of added stars and A^ r is the 
number of recovered stars. Fig. 4 shows the relation- 
ship between A and V magnitude. 
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Figure 5. CMD of the central region (r < 5') of NGC 
1245 with isochrones fitted by eyeball (red line) and \ 2 
minimum (green line). Black lines indicate the boundaries 
of main-sequence stars and giant box. 



Figure 7. Confidence regions of the joint variation in age, 
metallicity, distance modulus and reddening for NGC 1245. 
The three colors red, green, and blue represent the confi- 
dence limits of 2, 4, 6 <r, respectively, corresponding to 
Ax 2 = 4, 16, 36 



> 18 




Figure 6. CMD of the central region (r < 5') of NGC 
2506 with isochrones fitted by eyeball (red line) and \ 2 
minimum (green line). Black lines indicate the boundaries 
of main-sequence stars and giant box. 

3 PHYSICAL PARAMETERS 

We determined the basic physical parameters of the 
two open clusters by analyzing the color magnitude di- 
agrams (CMDs). We us ed Yonsei-Yale (Y 2 ) isochrones 
(|Demaraue et alj|20u4 ) to fit to the observed CMDs. 
We assumed E(V - I)/E(B -V) = 1. 25 and a total- 
to-se lective extinction ratio Rv of 3.2 |Cardelli et al.l 
1989). Initially, we adjusted the distance, age, metal- 
licity, and interstellar reddening by eyeball fitting 
(Fig. 5 and Fig. 6) but we determined the best fit 
parameters by applying x 2 minimization algorithm as 
described below. 




Figure 8. Confidence regions of the joint variation in age, 
metallicity, distance modulus and reddening for NGC 2506. 
The three colors red, green, and blue represent the confi- 
dence limits of 2, 4, 6 <r, respectively, corresponding to 
A X 2 = 4, 16, 36 

In order to reduce the effects of field star con- 
tamination, we constructed cluster CMDs using stars 
in the central region (R < 5') of the cluster. We also 
defined the boundaries of main-sequence stars, which 
are the blue and red edge along main-sequence stars. 
These boundaries are tightly selected in order to ex- 
clude binary sequence. According to evolution tracks, 
the blue edge of giant clump stars are more massive 
and more evolved than the red edge of the giant clump. 
And the giant branch of Y 2 isochrone represents the 



© 2012 RAS, MNRAS 000, [T|-?? 



6 Lee, Kang & Ann 



red edge of giant stars. So we defined a region of red 
edge of the giant clump. In Fig. 5 and Fig. 6, we 
showed the boundaries, and the red edge regions. The 
Y 2 isochrones provide an interpolation scheme to cal- 
culate isochro nes for an arb i trary age and metallicity. 

Following iBurke et al~l ((2004), for a given set of 
fitting parameters (age, metallicity, distance modulus 
and interstellar reddening), we defined the x 2 as 



(3) 



where the sum is over all stars selected for isochrone 
fitting and x? ls the contribution from star i which is 
defined as 



Wrr, 



V ]i + [(V - I) m -{V- /)„]? 



(4) 



Here, V and V — I) are magnitude and color of stars 
and the subscripts m and o stand for model and obser- 
vation, respectively. We assumed an equal weight for 
all stars used in the calculation of \ 2 m order to avoid 
overemphasizing the main-sequence turn-off where the 
systematic difference between theoretical isochrones 
are more pronoun ced (|Grocholski fc Saraiedinil [20031 ; 
IBurke et al.ll2004l ) although their photometric errors 
are the smallest. We performed a grid search over age, 
metallicity, distance modulus and reddening to find 
the best-fit parameters which minimize the \ 2 ■ I n ta- 
ble 3, we presents the ranges and intervals of param- 
eters used for NGC 1245 and NGC 2506. 

The confidence limits on the derived parame- 
ters are obtained b y scaling the resulting \ 2 statistic 
([Burke et al. 1120041 ) 



(5) 



where Xmin is the minimum x 2 and v is the number 
of degrees of freedom. In Fig. 7 and Fig. 8, we plot- 
ted the confidence regions for the joint variation in the 
isochrone parameters representing the A\ 2 = 4, 16, 36 
(red, green, blue), which correspond to 2, 4, and 6 
a error respectively. As shown in Fig. 7 and Fig. 8, 
the local minimum of A\ 2 is fairly well defined and 
the Ax 2 surface can be approximated by a paraboloid 
near minimum. In this approximation, the statistical 
errors (1 a errors) are the projection of the A\ 2 = 1-0 
extent of the paraboloid on the parameter axes. The 
best-fit parameters and their statistical errors (tJstat) 
are given in Table 4. As can be seen in Fig. 7 and 
Fig. 8, the statistical errors in the best-fit parame- 
ters are smaller in NGC 2506 than NGC 1245. The 
larger errors in NGC 1245 is caused by the broad 
main-sequence near the turn-off in NGC 1245. 

The statistical errors presented in Table 4 are very 
small. However, the uncertainties of the physical pa- 
rameters can be larger than these errors because they 
merely represent the statistical errors associated with 
the x 2 minimization. There are other sources of errors 
that affect the accuracy of the derived physical pa- 
rameters systematically. Possible causes of systematic 



errors are the uncertainties related to the photomet- 
ric calibration, assumed parameters such as the total- 
to-selective extinction ratio Rv, limiting magnitudes 
of observed data Vum and the theoretical models for 
isochrone fittings. In order to know the systematic er- 
rors in the derived physical parameters, we performed 
some experiments by varying the assumed parameters 
and theoretical models. We present the results in Ta- 
ble 4. 

From equation (1) and (2), systematic errors asso- 
ciated with the calibration of the photometry can be 
caused by uncertainties in the extinction and trans- 
formation coefficients and zero points. But we exam- 
ined only the effect of the uncertainties in the zero 
points (a z (V) and o z ) because the calibration coeffi- 
cients which are provided by CHFT web site are as- 
sumed to be well determined. We derived physical pa- 
rameters from the data with the zero point shifts of 
±Oz{V) = 0.05 and ±a z (V-I) = 0.10 and found that 
the systematic errors in E(B — V) and (V — My)o are 
relatively large. 

Since systematic uncertainties in Rv affect the 
cluster parameters, variations in Ry a lter the (V — I) 
color at fixed Av |Burke et al. I [20041 ). we performed 
X 2 minimization procedures to derive cluster parame- 
ters with Rv = 3.0 and Rv = 3.4 using the same grid 
search. We found that the systematic errors associ- 
ated with the Rv variation are similar to or less than 
the statistical errors, except for those in the distance 
modulus of NGC 1245. 

In Fig. 5 and Fig. 6, there is a non- negligible dis- 
crepancy between the observations and the theoretical 
isochrones for the lower main-sequence stars (V > 19). 
This effect is most serious for stars in 19 < V < 21 of 
NGC 2506. Due to this discrepancy, the derived clus- 
ter parameters can be changed if we assume different 
limiting magnitude for isochrone fittings. To estimate 
the uncertainties related to the discrepancy between 
the models and the observations, we performed the 
same procedures again to derive the cluster param- 
eters with limiting magnitudes of Vu m = 18.0 and 
Vu m = 21.0. The systematic errors related to the Vum 
are generally larger than the statistical errors and they 
are larger in NGC 1245. But, the systematic errors in 
(V-Mv)o are larger in NGC 2506. The reason for the 
larger uncertainties in the parameters of NGC 1245 is 
partly due to the broad distribution of stars near the 
turn-off and upper main-sequence of NGC 1245. The 
large systematic error in (V — Mv)o in NGC 2506 
seems to be caused by the maximum discrepancy be- 
tween the models and observations at 19 < V < 21 in 
NGC 2506. 

Since the cluster parameters derived from 
isochrone fittings significantly depend on the the- 
oretical models, w e used the Padova isochrones 
IIGirardi et ai] |2000D for comparison models. However, 
we did not apply x 2 minimization procedures because 
we do not have interpolation programs for the Padova 
isochrones. Rather, we applied eye-ball fittings for 
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Table 3. Isochronc grid parameters for NGC 1245 and 
NGC 2506 



Cluster 


NGC 1245 




NGC 2506 






Range 


Interval 


Range 


Interval 


age(Gyr) 


1.02 ~ 1.18 


0.005 


2.19 


~ 2.41 


0.005 


[Fe/H] 


-0.24 ~ 0.0 


0.005 


-0.29 


~ -0.16 


0.005 


(V - M v ) 


12.13 ~ 12.36 


0.005 


12.41 


~ 12.53 


0.002 


E(B - V) 


0.2 ~ 0.3 


0.005 


0.015 


- 0.070 


0.002 



Table 4. Physical parameters of NGC 1245 and NGC 2506 



Cluster 




E(B-V) 


(V - M v )o 


age(Gyr) 


[Fe/H] 


NGC 1245 


best v 2 . 


0.24 


12.25 


1.08 


-0.08 




zeropoint + tr z 


0.30 


12.20 


1.06 


-0.12 




zeropoint — u z 


0.21 


12.32 


1.08 


-0.09 




Rv = 3.0 


0.26 


12.33 


1.06 


-0.12 




R v = 3.4 


0.25 


12.18 


1.06 


-0.11 




V Hm = 18.0 


0.27 


12.33 


0.9 


-0.09 




V Hm = 21.0 


0.225 


12.23 


1.08 


0.00 




Padova (Eyeball) 


0.12 


12.45 


1.00 


0.07 




Y 2 (Eyeball) 


0.22 


12.35 


1.0 


0.0 




&stat 


±0.01 


±0.02 


±0.01 


±0.02 




O sys 


±0.05 


±0.12 


±0.09 


±0.06 




0~tot 


±0.05 


±0.12 


±0.09 


±0.06 




Previous 


0.26 2 , 0.28 s , 0.28 3 


13.2 2 ,12.0 3 ,12.4 5 


0.8 2 ,1.1 3 


-0.04 3 ,-0.05 6 






0.29 5 ,0.21 6 


12. 27 6 , 12.1 s 


0.9 5 ,1.05 6 


±0.06 s ,+0.14 9 


NGC 2506 


bef3t Xmin 


0.035 


12.47 


2.31 


-0.24 




zeropoint + a z 


0.08 


12.38 


2.30 


-0.24 




zeropoint — <j z 


0.00 


12.53 


2.32 


-0.24 




Rv = 3.0 


0.035 


12.47 


2.33 


-0.24 




R v = 3.4 


0.035 


12.46 


2.32 


-0.24 




V Hm = 18.0 


0.02 


12.33 


2.48 


-0.22 




V Hm = 21.0 


0.03 


12.32 


2.47 


-0.23 




Padova (Eyeball) 


0.02 


12.4 


2.5 


-0.32 




Y 2 (Eyeball) 


0.02 


12.4 


2.19 


-0.2 




Vstat 


±0.004 


±0.01 


±0.01 


±0.01 




O sys 


±0.04 


±0.08 


±0.16 


±0.06 




0~tot 


±0.04 


±0.08 


±0.16 


±0.06 




Previous 


0.1\0 ~ 0.07 4 
0.10 8 , 


11.7\ 12.6 4 ,12.95 7 
11.8 s 


1.6 ~ 2.0 4 ,1.8 7 


-0.32 ~ 0.07 4 ,-0.32 7 
-0.41 s 



MPurgathoferl dl964h. 2 ICarraro fc Pataj <1994l. 3 IWee fc Led lll996l), 4 lMarconi et al.l Jl997h . 5 ISubramaniaml 
l|2003l V. 6 iBurke et al. I l|2004l l. 7 lHenderson et aD J2007t > . s ljaned jl97Sft . 9 iLvngal {198% 



both the Y 2 isochrones and the Padova isochrones and 
compared the derived cluster parameters. As can be 
seen in Table 4, the systematic errors resulting from 
the theoretical models are much larger than the sta- 
tistical errors but comparable to the systematic errors 
associated with the assumed parameters. In general, 
the systematic errors associated with the theoretical 
models are larger in NCC 1245 than NGC 2506. The 
largest difference is observed in the age of NGC 2506. 

We calculated the systematic errors of cluster pa- 
rameters from the derived values in Table 4 by assum- 
ing that we explored typical values of Rv, Vu m and 
zero point shifts in the derivation of cluster parame- 
ters. Actually, we considered the standard deviations 



of the cluster parameters resulting from different sets 
of assumed parameters (zero points, Rv and Vu m ) and 
theoretical models as the systematic errors of the clus- 
ter parameters. Then, the total systematic error (a sys ) 
can be determined by a quadratic sum of the system- 
atic errors as 

o- S ys 2 = o~z 2 + o~(Rv) 2 + o(Vum) 2 + cr(iso) 2 (6) 

here a(Rv) and a(Vu m ) are systematic errors associ- 
ated with the parameters Rv and Vu m , respectively, 
and <j(iso) is the systematic error related to the the- 
oretical models. Then, the total error (o to t) can be 
written as 

o~tot 2 = o stat 2 + o- ays 2 (7) 
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where <T s tat is the statistical error related to the \ 2 
minimization. In Table 4, we present the total uncer- 
tainties in each cluster along with the statistical errors 
and the systematic errors. However, the actual uncer- 
tainties in the derived cluster parameters can be larger 
than these values since we did not consider the errors 
associated with the unresolved binaries and the field 
star contamination. Our est imates of the syste matic 
errors are similar to those of iBurke et al. I (|2004T ) who 
derived the cluster parameters by the same method 
as ours and the total uncert ainties are comparable to 
those of I Janes & Hod J 201 ll ) who derived the cluster 
parameters using synthetic CMDs. 



4 LUMINOSITY FUNCTION 

4.1 Construction of Luminosity Function 

Our observed area was wide enough to cover the field 
region as well as the cluster region. To derive the LFs 
of the open clusters, we first determined the cluster 
extent R c i by examining the surface density distribu- 
tion of the observed area. The radius at which the 
surface number density radial profiles of the stars be- 
comes constant was taken to be R c i . The field regions 
can then be defined as the region R > 2R c i. From 
this, we derived the field star population, which was 
required to correct the open cluster LFs. The cluster 
radii were R ct = 16' and 19' for NGC 1245 and NGC 
2506, respectively. We derived the LFs by counting 
the number of stars in each bin of magnitude Am = 1 
for the V-band photometry of the cluster region with 
the correction of field stars as follows: 



A 



Nj_Aci_ 
A Af 



(8) 



where N c i is the number of stars for each magnitude 
bin in the cluster region, Nf is that of the field region, 
and Ad and Af are the areas of the cluster and field re- 
gions, respectively. Here, we assumed that all the field 
regions have the same LF, but this is not usually true 
because there are fluctuations in the surface number 
densities of the field stars as well as inhomogeneous 
field star populations. 

In order to take into account the uncertainties 
arising from the fluctuations in different directions, we 
constructed LFs of the field regions surrounding the 
clusters (see Fig. 9 and Fig. 10 for NGC 1245 and NGC 
2506, respectively). The LFs of the NGC 1245 field 
regions were very similar to each other except for the 
large difference at Mv ~ 10. This large difference was 
caused partly by the different limiting magnitudes of 
the individually observed regions and partly by large 
photometric errors at the faintest magnitudes. While 
the four field LFs were nearly the same in NGC 1245, 
the field LFs in the south of NGC 2506 were much 
different from those in the north of NGC 2506. More 
specifically, the number of stars in the southern fields 
of NGC 2506 was larger than that in the northern 
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Figure 9. Luminosity functions of four field regions 
around NGC 1245. 




Figure 10. Luminosity functions of four field regions 
around NGC 2506. 



fields. Since NGC 2506 is located in the north of the 
Galactic plane, the difference seen in Fig. 10 was due 
to the density gradient along the Galactic latitude. 
Because the LFs of the open clusters were affected 
by those of the field regions, the dispersions in the 
field LFs are listed in Table 5 as a measure of the 
uncertainties caused by the correction. If we consider 
the errors of the cluster LFs as the dispersions in the 
field LFs scaled to the area of the cluster, they become 
larger than the Poisson errors, but this seems to be 
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Table 5. Dispersions in the LFs of field 
regions of NGC 1245 and NGC 2506 

Dispersion (stars/(arcmin) 2 ) 



M v 


NGC 1245 


NGC 2506 


-1.0 


0.007 


0.036 


0.0 


0.010 


0.046 


1.0 


0.004 


0.041 


2.0 


0.022 


0.090 


3.0 


0.033 


0.081 


1.0 


0.039 


0.117 


5.0 


0.024 


0.100 


6.0 


0.079 


0.193 


7.0 


0.137 


0.166 


8.0 


0.083 


0.128 


9.0 


0.142 


0.164 


10.0 


0.459 


0.135 


11.0 


0.192 


0.082 





Figure 12. Luminosity function of NGC 2506. Error bars 
represent the logarithmic dispersion of field LFs scaled to 
the cluster area. 
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Figure 11. Luminosity function of NGC 1245. Error bars 
represent the logarithmic dispersion of field LFs scaled to 
the cluster area. 
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more a mea ningful way to represe nt the errors in the 
cluster LFs (jPhelps fc Janedll993h . 

4.2 Characteristics of Luminosity Functions 

The LFs of NGC 1245 and NGC 2506 are shown in 
Fig. 11 and Fig. 12. The logarithmic error bar was 
scaled by area using Table 5 as follows: 

e+ = log(N + A x A cl ) ~ log(N) (9) 

e- = log(N) - log(N - A x Ad) (10) 

Here, e+ and e_ are the positive and negative errors, 
respectively, and A is the LF dispersion of the field 
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Figure 13. LFs compared with that of the solar neigh- 
bourhood. All LFs have been normalized to the value at 
My =4.0 



region as listed in Table 5. Our observing limit was 
Mv « 10, but near Mv « 10 for NGC 1245, the un- 
certainties increased dramatically because of the dif- 
ferent limiting magnitudes of each observing region. 

Stars brighter than My ~ 1 are thought to be 
turn-off and giant stars in NGC 1245 and NGC 2506. 
Comparing the bright parts of the LFs of NGC 1245 
and NGC 2506, the number of bright stars in NGC 
1245 is larger than that in NGC 2506. The reason for 
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the larger number of brighter stars in NGC 1245 is 
the younger age of NGC 1245. The age of NGC 2506 
is nearly twice that of NGC 1245. And the giant clump 
stars of Mv ~ are more developed in NGC 2506. 

The LF of NGC 1245 is nearly fiat. The flat LF of 
NGC 1245 could be caused by the dynamical evolution 
which makes the low-mass stars evaporate from the 
cluster because the LF of the cluster is expected to 
increase monotonically toward the faint magnitudes 
to at least M v « 10. The LF of NGC 2506 is similar 
to that of NGC 1245, but they do differ. While the 
LF of NGC 1245 is relatively flat between M v = 3 
and M v = 10, the LF of NGC 2506 shows a gradual 
increase toward fainter magnitudes. It seems unlikely 
that this difference is due to the errors in the field star 
correction because the dispersions in the field LFs are 
not large enough to erase the rising tendency toward 
the faint magnitudes. It seems worthy to note that the 
LF of NGC 2506 is derived from a slightly larger area 
than that of NGC 1245, which may include regions 
where evaporating stars are likely to be located, at 
least temporarily. The radial dependencies of the LFs 
of NGC 1245 and NGC 2506 show that the peak at 
M v = 10 in the LF of NGC 2506 is most pronounced 
in the LF derived from the outer regions. 

In Fig. 13, we compare the LFs of NGC 1245 
and NGC 2506 with that o f the solar neighbourhood 
(|Binnev fc Merrifieldl Il998h . Because turn off magni- 
tudes are about M v w 2.0 and 3.5 for NGC 1245 
and NGC 2506, respectively, we normalized the LFs at 
Mv = 4.0 to avoid evolutionary effects. In this figure, 
we can see a similar global pattern between the solar 
neighbourhood and the open clusters, but the number 
of bright stars in the LFs of NGC 1245 and NGC 2506 
is slightly larger than that of the solar neighbourhood, 
while the opposite is true for the number of faint stars. 
We believe this feature is caused by the evaporation of 
low-mass stars in NGC 1245 and NGC 2506 because 
they are old enough to expect this kind of dynamical 
evolution and mass segregation. 



5 MASS FUNCTION 

The mass function is a useful tool for studying the 
IMF and to understand the star formation process 
in general as well as the structure and evolution 
of a stellar system. In order to construct the mass 
function of NGC 1245 and NGC 2506, w e used the 
mass-luminosity relation of ISchaller et al.l lll992f) for 
stars with Mv < 7.0 and that of Henry fc McCarthy! 
(|l993h for M v = 7-17.6. Fig. 14 and Fig. 15 show 
the mass functions of NGC 1245 and NGC 2506 
derived by counting the stars in the mass interval 
Alog(M /Mo) = 0.1 inside R c i. The uncertainties in 
the mass function were derived from the dispersions 
of the field region mass functions, similar to the un- 
certainties of the LFs. We derived the slope of mass 
function 



NGC 1245 
Slope = -1.29 




0.2 o 

log(M/M ) 



Figure 14. Mass function of NGC 1245. Error bars rep- 
resent the logarithmic dispersion of field MFs scaled to the 
cluster area. The fitting range is log(M J M@) = —0.2 ~ 0.3. 




log(M/M ) 



Figure 15. Mass function of NGC 2506. Error bars rep- 
resent the logarithmic dispersion of field MFs scaled to the 
cluster area. The fitting range is log(M / Mq) = —0.2 ~ 0.1. 



dlogN{log{M/M Q j) 



(11) 



dlog(M/M Q ) 

by a least-squares fit, weighted by dispersions in the 
MFs of field regions. We obtained V = -1.35 ± 0.04 
for NGC 1245 and F = -1.26 ± 0.07 for NGC 2506 
with fitting ranges of log(M / 'Mq) — —0.2 ~ 0.3 and 
log(M/M Q ) = -0.2 ~ 0.1, respectively for NGC 1245 
and NGC 2506. 
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The derived F values of the two clusters were 
slightly shallower than th e slope of the solar neigh- 
bourhood IMF: T = -1.35 ljSalpeterlll955h . In the case 
of NGC 1 245, the present F v alue i s simi lar to those re- 
ported bv lCarraro fc Patatl (jl994) and ISubramaniaml 
(2003) who derived a similar but shallower slope than 
the Salpeter value. The slopes of the mass functions 
of NGC 1245 and NGC 2506 were thought not to be 
the same as the slope of the respective IMFs because 
of their age; evaporation of the low-mass stars results 
in a shallower slope. 



6 SUMMARY 

We conducted deep and wide VI photometry of old 
open clusters NGC 1245 and NGC 2506 based on CCD 
observations using the CFHT. The present photome- 
try is suitable for investigating dynamical evolution in 
open clusters because it is deep enough to reach down 
to My « 10 and wide enough to cover the entire re- 
gions of the two clusters including the field regions for 
an effective correction of the field star contamination. 

We derived the physical parameters of the two 
clusters using detailed isochrone fittings based on \ 2 
minimization: E(B-V) = 0.24 ±0.05 and 0.03 ±0.04, 
(V-Mv)o = 12.25±0.12 and 12.47±0.08, age(Gyr) = 
1.08 ±0.09 and 2.31 ±0.16, and [Fe/H]= -0.08 ±0.06 
and -0.24 ± 0.06 for NGC 1245 and NGC 2506, re- 
spectively. The quoted errors are the total errors which 
are quadratic sum of statistical errors associated with 
the Xmin an d systematic errors related to the assumed 
parameters and theoretical models. The present esti- 
mates of the physical parameters of two clusters are 
in good agreement with the previous estimates. 

We also derived the LFs of NGC 1245 and NGC 
2506. The LF of NGC 1245 shows a flat profile be- 
tween M v = 3 and My = 10, whereas the LF of NGC 
2506 displays a slight rise. This difference seems to be 
not due to the errors related to the field star correc- 
tion but due to the dynamical structures, which can 
differ because dynamical evolution can be different for 
clusters of a similar age as it is affected by the envi- 
ronment as well as by the internal properties of the 
cluster. We derived the mass functions of NGC 1245 
and NGC 2506, which show slightly shallower slopes 
than that of the solar neighbourhood IMF. The shal- 
lower slope is understandable if we consider that the 
derived mass function is the present day mass func- 
tion, which is thought to be different from the IMF 
of the clusters due to dynamical evolution, since the 
ages of NGC 1245 and NGC 2506 are old enough for 
mass segregation and evaporation of low-mass stars. 
We will discuss the dynamical properties, structures 
and halo of these two clusters in a forthcoming paper. 
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